Gravitino dark matter and baryon asymmetry from Q-ball decay in gauge mediation 
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We investigate the Q-ball decay in the gauge-mediated SUSY breaking. Q balls decay mainly 
into nucleons, and partially into gravitinos, while they are kinematically forbidden to decay into 
sparticles which would be cosmologically harmful. This is achieved by the Q-ball charge small 
enough to be unstable for the decay, and large enough to be protected kinematically from unwanted 
decay channel. We can then have right amounts of the baryon asymmetry and the dark matter of 
the universe, evading any astrophysical and cosmological observational constraints such as the big 
bang nucleosynthesis, which has not been treated properly in the literatures. 
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I. INTRODUCTION 

Production of dark matter and baryon number of the 
universe is very important in cosmology. In the gauge- 
mediated supersymmetry (SUSY) breaking, the gravitino 
is natural candidate of the dark matter, while the Afflick- 
Dine (AD) baryogenesis |l| is promising mechanism be- 
cause thermal leptogenesis confronts serious cosmological 
gravitino problem, for example. It is well known that the 
AD condensate transforms into nontopological solitons, 
Q balls, during the course of baryogenesis [2-6]. The 
charge of the Q ball in the gauge mediation could be large 
enough to be stable against the decay into nucleons, and 
the Q ball itself could be the dark matter of the universe 
|2J. In this case, however, nondetections in direct and/or 
indirect Q-ball searches tell us that it may be difficult to 
simultaneously obtain enough baryon asymmetry of the 
universe 1 . 

On the other hand, the small charge Q ball will de- 
cay into nucleons. At the same time, it could decay into 
gravitinos, which would have enough abundance to be the 
dark matter. This idea was considered in Refs.0, Q 2 . In 
the former, the authors considered relatively small charge 
Q balls which evaporate in the thermal bath, and imag- 
ined that the gravitino is directly produced from the Q 
balls. In this case, however, most of the charge decays 
(evaporates) at the temperature around the electroweak 
scale, and sparticles are created, finally leading to the 
production of the next-to- lightest SUSY particle (NLSP), 
which thermalizes immediately. The gravitino could be 
produced by the NLSP decay, but the abundance of the 
gravitino is too small to be the dark matter of the uni- 
verse, and the NLSP decay would destroy light elements 
during the big bang nucleosynthesis (BBN). 

In the latter reference, the formed Q ball is the gravity- 
mediation type, which is able to decay into NLSPs non- 
thermally. The NLSP decays into the gravitiono in the 



end, but they did not discuss the fact that the NLSP 
abundance is limited severely by the BBN constraints 
for m 3 /2 > 1 GcV, so that the amount of the produced 
gravitino should be much smaller than that to be the 
dark matter [9| 3 . 

In this article, we investigate a very simple scenario 
that the unstable Q balls decay mainly into nucleons, 
partially into gravitinos with small branching ratio, while 
kinematically forbidden to NLSPs almost throughout the 
decay process, in the gauge-mediated SUSY breaking. 
The crucial observation reveals that the right amounts of 
the baryon asymmetry and the dark matter of the uni- 
verse can be created through this process. The scenario 
is achieved by a couple of key ingredients. One is that the 
Q ball decays from its surface [llj so that the decay rate 
through the main channel has an upper bound, which is 
called the saturated rate, while the decay into gravitino 
is not saturated. Another is that a large enough Q ball 
forbids the decay into the NLSP, while allowing the de- 
cay into nucleons and gravitinos, having an appropriate 
lifetime to decay at the temperature of O(10 MeV). 

The whole scenario proceeds as follows. The AD field 
starts the oscillation (rotation) after inflation. The ro- 
tation has an oblate orbit in the potential. Due to the 
negative pressure, the AD field feels spatial instabilities, 
which grow into Q and anti-Q balls. The Q and anti-Q 
balls decay into nucleons and anti-nucleons to create the 
baryon asymmetry in the universe, while producing small 
amount of gravitinos which become the dark matter of 
the universe. NLSPs would be produced only at the very 
end of the decay process when the charge of the Q ball 
becomes small enough to open the kinematically allowed 
channel to the NLSP production[12|, and its abundance 
is small enough to evade any cosmological disasters such 
as spoiling the successful BBN. 

The structure of the article is as follows. In the next 
section, we show the features of the Q ball in the gauge- 
mediated SUSY breaking model, and the details of the 



1 More recent observational limits tighter than those applied in 
Ref. [3] may wipe out the allowed regions. 

2 Its possibility was mentioned in Ref. Q. 



3 Their model may work for a right-handed sneutrino LSP or if 
there is a small amount of R-parity violation llOH . 



decay process of the Q ball are described in Sec. IIII1 In 
Sec lIVl we estimate the abundances of the baryon num- 
ber, the dark matter, and the NLSP. We impose BBN 
constraints on the NLSP abundances to derive the al- 
lowed gravitino mass range in Sec|V] In Sec lVIi we show 
the region in the model parameter space where the suc- 
cessful scenario exists. Section IVIII is devoted to our con- 
clusions. 



II. Q BALL IN GAUGE MEDIATION 

The AD field consists of some combination of squarks 
and sleptons. The potential is flat in the SUSY-exact 
limit, but due to the breaking of SUSY in gauge media- 
tion, it is lifted as[13| 
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where m^ ~ O(TeV) is a soft breaking mass, and M F and 
Ms are related respectively to the F and A components 
of a gauge-singlet chiral multiplet S in the messenger 
sector as 



the field value at the center of the Q ball are related to 
the charge Q as 
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respectively. 



III. Q-BALL DECAY 

A Q-ball decay takes place if some decay particles carry 
the same kind of the charge of the Q ball, and the mass 
of all the decay particle is less than Mq/Q. In this case, 
the Q-ball decay rate Tq has an upper bound, which we 



call the saturated rate Tq 



Tr> < r (sat) 



1 



Q 



J Q 



Q 



Q 192tt : 



■ATrRi 



24-y/2 



M F Q 



-5/4 



(10) 



Mf. 



(4tt) 



-(^5> 2 , 



M S = (S). 



(2) 



Here g generically stands for the standard model gauge 
coupling, and Mp is allowed in the following range: 



where we use Eqs.® and © in the last equality. The de- 
cay rate saturates typically when f e s4>Q ^ wq , where the 
elementary process has an interaction of £j nt = f e s<fyipil). 
On the other hand, for f c s4>Q <C uiq, the decay rate is 
given bypjj 
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where m 3 / 2 is the gravitino mass and Mp 
GeV is the reduced Planck mass. 

The AD field starts the oscillation when H ~ M F /(j> osc , 
where OSC is the field amplitude at the onset of the oscil- 
lations. The field fluctuations grow exponentially due to 
the negative pressure during the helical motion, and the 
field transforms into Q balls. The charge of the formed 
Q ball is estimated asp 
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where (3 ~ 6 x 1CP 4 for a circular orbit (e = 1), while j3 ~ 
6 x 10 -5 for an oblate case (e < 0.1). Here e represents 
the ellipticity of the field orbit. The charge Q is just the 
^-numbers, and relates to the baryon number of the Q 
ball as 



B = bQ, 
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where b is the value of the baryon number carried by a 
* particle. For example, b — 1/3 for the udd direction. 
The mass, the size, the rotation speed of the field, and 
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Since we are interested in the case that the Q ball can 
decay into nucleons (Mq/B > m^, where mjv — 1 GeV 
is the nucleon mass) and, at the same time, the decay into 
NLSPs is kinematically prohibited (Mq/Q < mNLSp), 
the main decay channel is the decay into nucleons. We 
thus consider the Q-ball charge in the range Q cl - < Q < 
Qd, where 
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Notice that the NLSP is produced at the latest moment 
of the decay, only after the charge of the Q ball becomes 
less than Q cr during the course of the decay [12l|. 

The elementary process of the main channel is squark 
+ squark — > quark + quark via (heavy) gluino exchange, 
whose rate is estimated as|14] 
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where m~ g is the gluino mass, and £ ~ \Vc 
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is a pos- 



sible CKM suppression factor (10 3 < |Vckm| ^ 1 
Thus, we obtain the effective coupling as 
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the main decay channel is saturated so that we must use 
the saturated rate Tq ' for the decay. Thus, the Q ball 
decays when the cosmic time becomes the lifetime of the 
Q ball: t ~ 1/Iq . The Q ball decays when the universe 
is radiation-dominated, and the cosmic temperature at 
the decay is estimated as 
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where N* is the degrees of freedom at the corresponding 
temperature, and set to be 10.75 here. 

On the other hand, the Q-ball decay into gravitinos is 
not saturated. The decay rate of the elementary process 
squark — > quark + gravitino is given by 
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we have 
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where we use Eqs.® and ((9|) in the first equality, and 
Eq.(|3j) in the second equality. Therefore, the branching 
ratio of the decay into the gravitino is estimated as 
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IV. BARYON NUMBER, GRAVITINO DARK 
MATTER, AND NLSP ABUNDANCES 

The number densities of the baryon number, the grav- 
itino dark matter, and the NLSP are respectively related 



to the number density of the AD field as 
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WMAP seven-year data tells us pvM/pb = 4.94 + °;g3^5|. 
Since the ratio of the gravitiono dark matter and the 
baryon densities is 
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Therefore, the orbit of the AD field should be oblate, 
which is natural in the gauge-mediated SUSY breaking 
models. As mentioned above, this leads to the produc- 
tion of both Q and anti-Q balls to cancel out the baryon 
number to be small enough, compared to the gravitino 
dark matter. 

Baryon number is created when the AD field starts the 
rotation. It can be estimated as 
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are used in the last equality in the second line, and 
0OSC = /3~ 1/4 M F Q^ 1 / 4 from Eq.® is used in the last 
line. Inserting Eq. (|2l))) . we obtain 
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This provides the relation among Q, Mp, and Trh- 
On the other hand, the NLSP abundance is estimated 

as 
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where Y 3 / 2 = n 3 / 2 /s, and we use P3/2/ ' Pb — 5 in the 
second equality, and EdS-PS]). (PHI) , (fltfl) . and d2U) in 
the third one. We can now consider the limits of the 
abundance from BBN constraints in the next section. 



V. BBN CONSTRAINTS 

The abundance of the NLSP is limited by BBN con- 
straints, since the decay of the NLSP affects the abun- 
dances of light elements. It depends on the gravitino 
mass and the species of the NLSP how it has an influence 
on them. Let us first estimate the expected abundance 
of the NLSP in the limiting case, which means that the 
Q ball decays just in time for BBN (T D =1 or 3 MeV), 
and the thermally produced gravitinos do not overdom- 
inate. The latter condition is rephrased in terms of the 
constraint on the reheating temperature as[16j 
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ing Eqs. fTT]) and (f^5|) . we can eliminate Mp and Q from 
Eq.(f2"TJ|). which results in 
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where Eq. (f30"|) is used in the last inequality. In Fig|T] it 
is shown in red lines, above which is allowed. Also shown 
are the BBN constraints for bino (blue), stau (pink), and 
sneutrino (green) NLSPs, where upper right regions are 
excluded Q. In addition, we plot the dotted line as the 
NLSP abundance that the gravitinos produced by the 
NLSP decays has a right amount of dark matter: 
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Therefore, the NLSP abundance should be larger than 
the limit Eq.([3"Tj). and smaller than either excluded re- 
gions by BBN constraints in each NLSP species, or the 
upper limit that the produced gravitinos from the NLSP 
decay do not overdose the universe [Eq. (|3"2"j) ]. One can 
see that the gravitino mass above GeV is excluded. The 
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FIG. 1: NLSP abundance Eq.(|3T) for T D = 1 (3) MeV, de- 
noted as the lower (upper) red line, for (a) ?tjnlsp = 300 GeV 
and (b) winlsp = 100 GeV. BBN bounds are shown for bino 
(blue), stau (pink), and sneutrino (green, only in (a)), taken 
from Ref. 9] . Dotted line represents the upper bound of the 
abundance of NSLP that decays into gravitinos, above which 
overdose the universe [Eq. (|32[) ]. The allowed abundance is in 
between the red line and BBN bounds or the dotted line. 



allowed region is typically m 3 /2 ^ GeV for tonlsp = 300 
GeV, and m 3/2 < 10~ 2 GeV for to N lsp = 100 GeV. The 
detail will be found in Table |U 



VI. SUCCESSFUL SCENARIO 

As shown in the previous section, the allowed param- 
eter space is restricted by BBN constraints severely, but 
there is still a good chance to have successful scenario. 
Here we estimate the allowed region in Q-ball parame- 
ters, namely in Q — Mp plane. 

Equation (|2"8|) tells us the iso-TRH lines in Q — Mp 
plane, which reads 
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TABLE I: Upper limit of m 3/2 in GeV. 







T D = 1 MeV 


T D = 3 MeV 


300 GeV 


bino 


0.48 


3.9 x 10" 2 


300 GeV 


stau 


1.6 


0.86 


300 GeV 


sneutrino 


0.91 


0.66 


100 GeV 


bino 


2.5 x 10" 3 


1.9 x 10" 3 


100 GeV 


stau 


5.3 x 10" 2 


2.1 x 10 -2 



In the same manner, we can obtain iso-NLSP density 
lines from Eq. (|29l) as 
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Upper limit of the reheating temperature depends on 
the gravitino mass 7713/2 in such a way as in Eq. (|30[) . and 
once we fixed m 3 / 2 we can restrict the NLSP abundances 
from FigfT] Therefore, for the fixed 1113/2, the allowed 
region is above the line (|33p with the largest possible 
Trh, and between the lines (|34p with the lower and upper 
NLSP abundances. In addition, the Q ball must decay 
before the BBN time so that the temperature at the decay 
Td should be larger than 1 MeV (or 3 MeV for more 
conservative limit), represented as line (|3"5)) . 

We plot these lines for ttinlsp — 300 GeV in Q — Mp 
planes in FigHfor m 3/2 = 1 GeV, 0.1 GeV, 0.01 GeV, 
and 1 MeV. Red lines denote the T-q contours, light blue 
lines are the Xrh contours, and dark green lines represent 
the NLSP abundance contours. Also shown in blue are 
the mass per charge Mq/Q lines: 
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above which is prohibited for the Q ball to decay into 
those particles with mass larger than Mq/Q. In par- 
ticular, the Q-ball decay with baryon number is stable 
above the upper blue line denoted by Mq/(b Q) — 1 GeV 
[equivalent to Eq.([l3"|)]. and the 300 GeV NLSP is only 
produced once the charge of the Q ball becomes below 
the lower line denoted as Mq/Q = 300 GeV [equivalent 
toEq.(H2D]. 

Only the stau NLSP is allowed in the region shown in 
faint red in FigfSJa). In Fig[2tb), there are three territo- 
ries, where the faint (T D > 1 MeV) and the dark (T D > 3 
MeV) red are applicable to the stau and the sneutrino 
NLSPs, while the blue is for the bino NLSP. On the other 



hand, all three species (bino, stau, and sneutrino NLSPs) 
are relevant for the allowed region in Fig^c), which cor- 
responds to the NLSP abundance between the red lines 
[Eq.dSJ)] and the dotted line [Eq.fl35J)] in FigOJa). We 
can see in Fig.^d) where 77i 3 / 2 = 1 MeV that the allowed 
region disappear completely for Td > 3 MeV, and only 
a tiny area remains for Td > 1 MeV. Thus, we obtain 
the lower bound of the gravitino mass that could result 
in successful scenario as m 3 / 2 > 1.1 X 10~ 3 (6.9 x 10~ 4 ) 
GeV for T D > 3 (1) MeV. 

Notice that the allowed regions reside in such parame- 
ter space that the Q-ball charge is by far large enough so 
that the Q ball survives from the evaporation in thermal 
bath, where the Q ball evaporates if the charge is less 
than d 
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which is shown by green dashed lines in Figf2] On the 
other hand, black dotted lines in these figures denote the 
boundary above which the Q ball is the gauge-mediation 
type: 
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Also notice that those Q balls in the allowed region in 
Fig. [2ja) (7713/2 = 1 GeV) will form when the field starts 
rotation at the field amplitude larger than ~ Mj./m 3 / 2 , 
where the gravity-mediation effects dominate over the 
gauge-mediation ones. Depending on the one-loop_po- 
tential, the formed Q ball could be the delayed type [y| or 
the new type [18| which later transforms into the gauge- 
mediation type, or the gauge-mediation type may form 
directly. 



VII. CONCLUSIONS 

We have investigated the scenario of simultaneous pro- 
duction of the baryon asymmetry and the dark matter 
of the universe through the Q-ball decay in the gauge- 
mediated SUSY breaking model. This is simply achieved 
by the Q ball with charge small enough to decay into nu- 
cleons to create baryon number of the universe. We have 
calculated the branching ratio of the gravitino to find 
that it is small. Therefore, in order to obtain the ob- 
served baryon to dark matter abundance ratio, the AD 
field should rotate in the oblate orbit so that both Q and 
anti-Q balls are created to suppress the baryon number 
compared to the gravitino dark matter. 

At the same time, the charge should be large enough to 
kinematically forbid the decay into NLSPs. This evades 
the serious BBN constraints on the NLSP abundance, 
which was not properly taken into account in the litera- 
tures. NLSPs are actually produced at the very end of 
the Q-ball decay when the charge becomes small enough 
to open the corresponding decay channel. 



Including all these considerations, we have used the 
BBN constraints for the bino, stau, and sneutrino NLSPs 
in Ref.0), to get the allowed range of the gravitiono mass. 
We have found that typically m 3 / 2 ^ 1 (10~ 2 ) GeV for 
"inlsp = 300 (100) GeV. 

We have also shown the allowed region in Q — Mf 
plane for ttinlsp = 300 GeV. We thus have found the 
production of right amounts of both the baryon asym- 
metry and the gravitino dark matter of the universe is 
successful when the Q ball has the charge Q ~ 10 23 for 
M F ~ 10 6 - 10 7 GeV and m 3/2 ~ 10~ 2 - 10" 1 GeV. It is 
naturally achieved, for example, for the n — 6 udd direc- 
tion. In this case, the nonrenormalizable superpotential 
W^nr = (udd) 2 /Mp determines the field amplitude at the 
onset of the oscillation <^> osc ~ (MpMp) 1 / 5 , and also leads 



to the A-teim of the form Va ~ am 3 / 2 M / NR + h.c, and 
the scenario is successful for Mp ~ 10 7 GeV, m 3 / 2 ~ 0.01 
GeV, and a ~ 0.1, which leads to (f> osc - 7 x 10 13 GeV, 
and Q ~ 10 23 . 
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FIG. 2: Allowed region in Q — Mf plane for (a) m 3 / 2 = 1 GeV, (b) m 3 / 2 =0.1 GeV, (c) m 3 / 2 = 0.01 GeV, and (d) m 3 / 2 = 1 
MeV. In all figures, red lines denote iso-TD contours from 1 MeV to 1 GeV from the top to the bottom, respectively, and blue 
lines denote the limit of the mass of the particle which the Q ball can decay into: the upper is for nucleons and the lower is for 
the 300 GeV NLSP. Here we take b = 1/3. Gray hatched region is out of the range of Mf [Eq.©]. 

(a) Only the stau NLSP has allowed region. Light blue lines represent the Trh contours for Trh = 3 x 10 7 GeV (the lower) 
and 3.5 x 10 5 GeV (the upper). The upper and lower dark green lines indicate Pnlsp/s = 2.3 x 10~ 9 GeV and 1.0 x 10 -8 GeV, 
respectively. 

(b) The faint and the dark red regions are applicable to the stau and the sneutrino NLSPs for To > 1 MeV and To > 3 MeV, 
respectively, while the bino NLSP is only allowed in the blue region. Light blue lines represent the Trh contours for Trh = 3x 10 6 
GeV (the lower), 74 GeV (the middle), and 2.7 GeV (the upper). The dark green lines indicate pnlsp/s = 4.9 x 10~ 9 GeV, 
8.8 x 10~ 9 GeV, 4.4 x 10~ 8 GeV, and 1.5 x 10~ 6 GeV, from the top to the bottom, respectively. 

(c) The faint and the dark red regions are applicable to all the NLSP species for To > 1 MeV and To > 3 MeV, respectively. 
Light blue lines represent the Trh contours for Trh = 3 x 10 GeV (the lower), 64 GeV (the middle), and 17 GeV (the upper). 
The dark green lines indicate Pnlsp/s = 1.0 x 10~ 8 GeV, 9.4 x 10~ 8 GeV, and 1.5 x 10" 5 GeV, from the top to the bottom, 
respectively. 

(d) Allowed region is appicable to all the NLSP species. Light blue lines represent the Trh contours for Trh = 3 x 10 4 GeV (the 
lower) and 1 x 10 4 GeV (the upper). The upper and lower dark green lines indicate Pnlsp/s = 2.3 x 10~ 8 GeV and 2.0 x 10~ 7 
GeV, respectively. 

In all figures, green dashed line shows the limit of the charge evaporated in thermal bath, and Q ball is the gauge-mediation 
type above the black dotted line. 



